Coherence properties of spinor condensates at finite temperatures 
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We consider a spinor condensate of *^Rb atoms in its F = 1 hyperfine state at finite temperatures. 
Putting initially all atoms in mp — component we find that the system evolves into the state 
of thermal equilibrium. This state is approached in a step-like process and when established it 
manifests itself in distinguishable ways. The atoms in states mF = +1 and thf = — 1 start to rotate 
in opposite directions breaking the chiral symmetry and showing highly regular spin textures. Also 
the coherence properties of the system changes dramatically. Depending on the strength of spin- 
changing collisions the system first enters the stage where the mp = -1-1 and rriF = — 1 spinor 
condensate components periodically loose and recover their mutual coherence whereas their thermal 
counterparts get completely dephased. For stronger spin changing collisions the system enters the 
regime where also the strong coherence between other components is built up. 



Since a ^■^Na atoms Bose-Einstein condensate was cre- 
ated in a pure optical trap [1;] it has become possible 
to study the spinor properties of alkali dilute quantum 
gases. In such (spinor) condensates the atoms can change 
their spins upon collisions and still remain trapped. 
Therefore, in spinor condensates the population of their 
components is allowed to change in time. This new de- 
gree of freedom distinguishes the spinor condensates from 
mixtures of quantum gases. 

There has been recently a growing interest in such sys- 
tems. Different ground-state phases (ferromagnetic or 
polar) have been observed experimentally in = 1 con- 
densates [3, 0]- In particular, it has been shown that 
the ground state of rubidium condensate exhibits a fer- 
romagnetic behavior whereas the sodium system reveals 
the polar phase. For a F — 2 manifold the set of possible 
ground states of spinor condensates is even richer, includ- 
ing the so-called cyclic phase [5] . Experiment proved that 
under the natural conditions a ^^Rb condensate shows a 
polar behavior Recently, a first F = 3 (^^Cr) con- 
densate has been achieved [6| and many new spin phases 
were predicted for this system Q . 

There has been also recently an interest in the dynam- 
ics of the spinor condensates, including the coherent os- 
cillations between various components formation of 
domains Q, magnetically tuned resonances 0], and pre- 
dicted theoretically the effect of breaking of chiral sym- 
metry Some theoretical effort has been made to 
explore the dipolar interactions in addition to the short- 
range (van der Waals) interactions. In Ref. 0, [HI the 
coupling between the spin and orbital motion, inherently 
built in the dipolar interaction, is investigated and the 
phenomenon of the Einstein-de Haas effect is exhibited. 
As has been shown in Ref. II4I this effect might be ob- 



servable even in a very weak dipolar systems such as ru- 
bidium condensate provided the experiment is performed 
under resonance condition. 

In this Letter we investigate temperature properties of 
a spinor condensate of ®^Rb atoms in F = 1 hyperfine 
state. The influence of finite temperature on spinor evo- 



lution has been studied recently experimentally [13| and 
theoretically although in a quasi-one-dimensional system 
[l3 |. In the second quantization notation, the Hamilto- 
nian of the system is given by [l^ 



H = I d^r 



Y^l(i")V'l(r)V'i(r)V'i(i-) 



+ yV'i(r)it(r)F,,Ffc,V',(r)^Kr) 



(1) 



where repeated indices have to be summed over the val- 
ues mp = +1,0,-1. The field operator ^l^i{r) annihi- 
lates an atom in the hyperfine state \F = l^mp = i > 
at point r. The Hamiltonian consists of the kinetic en- 
ergy and the trapping potential (iJo — ~jm'^'^ + ^tr, 
where M is the mass of an atom) and two terms which 
describe the spin-independent and spin-dependent parts 
of the contact interactions, respectively. Coefficients co 
and C2 can be expressed with the help of the scatter- 
ing lengths flo and 02 which determine the collision of 
atoms in a channel of total spin and 2. The appropri- 
ate formulas are given by cn = 47r?i^(ao -I- 2a2)/3Af and 
C2 = 47rft2(a2 - ao)/3M According to Ref. [3, 

ao = 5.387 nm and 02 = 5.313 nm. Moreover, F are the 
spin-1 matrices. 

Within the classical fields approximation [l3] one re- 
places the field operators ipi (r) by the classical wavefunc- 
tions ipiii^)- The equation of motion for these wavefunc- 
tions is written as 



(2) 



The diagonal part of Ti is given by Hn = Hq + (cq + 
C2) \ipi\^ + (co + C2) IV^oP + (co - C2) Tioo = Ho + 

(co + C2)|V'lP + Co|V'Op + (co + C2)|V'-lP, H-1-1^Hq + 

(co - C2) IV'iP + (co + C2) l^oP + (co + C2) iT/i-ip. The 
off-diagonal terms describe the collisions which do not 
preserve the projection of spin of each atom (however, 
the total spin projection is conserved) and equal Hio — 
C2 ■0-1^0, '^0-1 = 02^1^1- Moreover, Hi-i = 0. 
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To study the properties of a spinor condensate at finite 
temperatures we prepare initially an atomic sample in the 
niF — component Q ^^'^ allow the system to evolve. 
The Bose gas is confined in a pancake-shaped trap with 
the aspect ratio (3 — LUz/tOr — 80 and the radial frequency 
ciJr = Svr X 100 Hz. To simplify calculations, we solve nu- 
merically the Eq. 12]) on a two-dimensional grid using 
rescaled scattering lengths (i.e., multiplied by a factor 
V^/3/27r). Obviously, spin changing collisions allow tran- 
sitions to other Zeeman states {mp — ±1) and hence 
the interplay between the spin dynamics and the ther- 
mahzation processes can be studied. In Fig. [1^ we plot 
the kinetic energy of each spinor component as a func- 
tion of time for the initial thermal noise introduced to 
the mp = component resulting in its ~ 20% depletion. 
We find that after some time the kinetic energy becomes 
equally distributed among spinor components (see Fig. 
H]), suggesting the system reached an equilibrium. The 
time needed for that depends on the initial amount of 
the thermal noise (the higher thermal energy the longer 
relaxation time). However, approaching this equilibrium 
is rather a step-like process not a uniform one. This 
is because the atomic cloud becomes fragmented and the 
overlap of different spin components changes significantly 
in time (Fig. [T^, green curve). The steps in Fig. [Iji seem 
to correspond to the maximal overlap. 

Spin changing collisions between atoms in mp — 
state produce equally atoms in states mp = ±1. Initially, 
both condensed and thermal atoms go to the mp ~ ±1 
components and back (see Fig. [T]d). This behavior 
changes after the equilibrium is established, i.e., after 
0.6 s. From now on, the fraction of thermal atoms in each 
component becomes 1/3 (Fig. [T]3, main frame). Since 
the total number of thermal atoms is also constant the 
system reaches the thermal equilibrium meaning all spin 
components in the thermal cloud are equally populated 
(already observed experimentally [HI). At the same time 
a non-trivial spin dynamics for condensed atoms is ob- 
served (Fig. [T]3, inset). Thus, the spinor condensate 
decouples from the thermal cloud. Also, since the con- 
densed atoms do not populate equally spinor components 
(as the thermal atoms do) it means that the single mode 
approximation is inappropriate at finite temperatures. 

To discuss the coherence properties of spinor conden- 
sates it is necessary to have in mind the stability dia- 
gram for such a system. This diagram is formed as a 
result of the Bogoliubov-de Gennes analysis of Eq. ([2]) 
in two-dimensional space and is considered with respect 
to the parameters (c2,co). More precisely, the stabil- 
ity of mp = component against the excitations in 
mp = ±1 states is investigated with particular empha- 
sis on modes of types (x e™"^, where m = ±1,0 (see 
Ref. 10] ). This analysis proves the existence of the re- 
gions where all m — ±1, modes are stable, only modes 
m = ±1 are unstable, and all m — ±1,0 modes are 
unstable. The coherence properties strongly depend on 
which region the systems parameters belong to. For ex- 
ample, if modes m = ±1 get unstable the superposition 
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FIG. 1: (color online), (a) Kinetic energy of spinor com- 
ponents as a function of time for mp = (red color) and 
mp ~ ±1 (black color) Zeeman states. The green line is 
the overlap between the +1 and —1 components defined by 
the scalar product |(^/j+i, i/)_i)|. (b) Fraction of total ther- 
mal (main frame) and total condensed (inset) atoms in each 
spinor component (red for and black for ±1 states). 



of (X e"^ and cx e^"^ states is excited in mp = 1 (as 
well as in mi? = —1) component. Assuming both modes 
equally contribute to this superposition one should ex- 
pect the fragmented (and consisted of two parts) densities 
in mi? = ±1 components. Examples are given in Fig. [21 
This fragmentation and circulation of atomic clouds are 
responsible for the coherence properties discussed below. 
Note also that the thermal atoms cause the granularity 
of densities since decreasing the level of initial thermal 
noise leads to very smooth densities. 
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FIG. 2: (color online). Density of -1-1 (left frame) and — 1 
(right frame) components for C2 = —1.0 x 10"'^, co = 0.22, 
and Ekin{0) = 9 x 10^ osc. units at 0.8 s. 

To investigate in detail the coherence properties of 
spinor condensates we calculate the one-particle density 
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matrix by averaging the product \E'*(r, t)4'(r', t) (with 
the spinor classical field ^(r,p being the solution of Eq. 
(121)) over the finite volume Next, this matrix is diag- 
onalized and an eigenvector corresponding to the domi- 
nant eigenvalue represents the spinor condensate wave 
fmiction. Let us denote it by ipo{r,t). Then the thermal 
component of the spinor classical field can be defined as 
-0"^ = 5* — {(po, 5*) (fio (here, the vector ipo is normalized 
to unity). Thermal part of the classical field lives in a 
space orthogonal to the condensate wave function. Now 
we introduce the spin density matrices for the condensate 
and the thermal part as below 



P^J 



(3) 



These matrices tell us about the inherent coherence of the 
spinor condensate and the thermal cloud. It is convenient 
to consider the normalized spin density matrices pij = 
Pij I y/PiiPjj- The diagonal elements of such matrices are 
equal to one and their off-diagonal elements determine 
the degree of coherence between various components of 
spinor condensate or the thermal cloud. When the off- 
diagonal element takes the value of the order of one (note 
that the maximal value of | pij \ equals 1 ) the coherence is 
strong whereas opposite means that the components are 
completely dephased. 

Coherence properties of '^^Rb spinor condensate are 
displayed in Fig. [3l where we show time dependence of 
spin density matrix elements corresponding to different 
regimes of parameters according to the stability analysis. 
When the coupling constant C2 is negative (ferromagnetic 
phase) and large enough to develop the rotational in- 
stability (i.e., when modes m = ±1 become unstable) 
the strong coherence between the -1-1 and —1 spinor con- 
densate components builds up (upper frame). This co- 
herence is maximal while the system is approaching the 
thermal equilibrium and remains strong afterwards, how- 
ever, in a way that it is periodically lost and restored. 
Surprisingly, even thermal cloud in components and 
— 1 is fully coherent during the way of the system to its 
equilibrium but later on (after 0.6 s) the thermal compo- 
nents get completely dephased. In this particular regime 
of parameters (c2,Co) we observe only a development of 
coherence between +1 and —1 components. This can be 
understood by comparing the overlaps |('0-i-i, cind 
1(0+1, ■f/'o) I during the evolution - the latter one turns out 
to be neghgible. 

However, the situation changes when the strength of 
spin-changing collisions C2 becomes large enough (and 
still negative) to develop an instability leading to signifi- 
cant population oi mp = ±1 Zeeman states with non- 
rotating atoms (i.e., when the mode m = becomes 
unstable). In the case of Fig. the system reaches 
its equilibrium already after 0.2 s and then the overlap 
I (■0+1, 00)1 gets even larger than the KV^+i, ■0-i)| one and 
it results in a development of coherence between -1-1 and 
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FIG. 3: (color online). Spin density off-diagonal elements 
for (a) (-1-1,-1) components (c2 = -1.0 x 10"^, Co = 0.22, 
and Eki„{Q) = 32 x 10^) and (b) (+1,0) components (c2 = 
-2.29 X 10"^ Co = 0.50, and Ekir^[Q) = 22 x 10^) for the 
condensate (black line) and the thermal cloud (red line) as 
a function of time. The overlaps |(i/)+i, i/>_i)| and K^/j+i, ?/>o)| 
are plotted in green in (a) and (b), respectively. 



components. The same is true for the pair of —1 and 
states. Therefore, in this range of (c2,co) parameters all 
components of spinor condensates remain mutually co- 
herent. Note that the thermal parts in +1 and states 
are always dephased. 

Both frames show also the overlaps: |(-0+i, V^-i)] in 
(a) and |('0+i,0o)| in (b) case in green color. In the 
equilibrium the spin density off-diagonal elements follow 
the corresponding overlaps. It clearly indicates that the 
coherence is built as a result of many atomic collisions. 
They allow to establish the relative phase of the compo- 
nents. In case (a) in the equilibrium the phase of the off- 
diagonal element pi_i becomes fixed within the intervals 
which correspond to the maximal overlap but changes 
going from one maximum to the other. On the other 
hand, in (b) the phase of pio gets constant just after the 
system reaches the thermal equilibrium (i.e., after 0.2s). 

Another peculiar feature of spinor condensates re- 
vealed after the system reaches its equilibrium is breaking 
of the chiral symmetry (at zero temperature this phe- 
nomenon was already reported in 10]). Since the modes 
m — ±1 get unstable the atoms in components -1-1 and 
— 1 start to rotate in opposite directions (the total or- 
bital angular momentum projection of colliding atoms is 
preserved by the contact interactions) on the onset of the 
equilibrium. No spontaneous rotation is observed when 
the system is approaching the equilibrium state. Fig. 0^ 
(curves black, red, and green) shows the typical behav- 
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ior of the angular momentum of all components when 
co/c2 — —216.4 (this is the natural ratio for F = 1 *^Rb 
atoms). In components +1 and —1 the atoms rotate 
alternatively clockwise and counterclockwise. This prop- 
erty changes when co/c2 differs strongly from —216.4. In 
this case components rotate only in one direction perma- 
nently breaking the chiral symmetry (blue curve in Fig. 
2^). Note, that unlike in [10] this permanent rotation is 
obtained here in the conservative, hamiltonian dynamics. 



(a) 



(b) 




0.4 0.6 
TIME (s) 



7.5 



-7.5 



■ ■ ■ ~ ^ M=50.-54 s 

. . ♦ * V *N \ * . . 

, , ^^ \ f f , 

■ ' i \,\ 

' I i \ \ >» - " ■ 

■ ' 1 ^ N ■• - ■ 
■ . \ \ >t * - ■ 



■ - " * - t=e.74 

f ^W^^- 

' ' ^* w ' * 

\ t 

■ ■ - s \ I. > 



• • - ^ t=e.97 s 

V -N V X « . 



-7.5 



7.5 



FIG. 4: (color online), (a) Angular momentum of components 
rriF = +1 (black curve), mp = ^ (red curve), and mp ~ ~1 
(green curve) as a function of time for parameters C2 = —1.0 x 
10"^, CO = 0.22, and Ekin{0) = 9 x 10^. Blue curve shows 
the angular momentum of ttif = +1 state in the case when 



co/c2 / -216.4 (C2 



-4.06 X 10" 



Co 



0.25). (b) Spin 



textures at particular instants of time for parameters (c2,co) 
as in the first case of part (a). The length of the vector (as well 
as the color) represents the quantity {< >^ + < Fy > 
and the size of each frame is 15 x 15 oscillatory units. 
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It turns out that the periods the atoms in mp = +1 or 
mj7 = — 1 states rotate in one direction (black and green 
lines in Fig. [4^) coincide with those when the off-diagonal 



spin density matrix element |pi-i| takes large values, 
i.e., when the mutual coherence between tof = +1 and 
tof = — 1 components is built up. In these time inter- 
vals the system develops also highly regular spin textures. 
Examples are given in lower panel in Fig. [4] where we 
plot the vector field composed of the average values of 
the total spin at particular instants of time. The first 
(from the left) frame shows the spiral galaxy-like pattern 
which persists within the interval (0.52 — 0.62) s. Later 
on a quadrupole-like texture is developed and lasts from 
0.74s to 0.81s (the second frame). The third plot shows 
domain-like pattern. On the other hand, for the case 
when atoms rotate in one direction (blue curve) the sys- 
tem exhibits all the time the same structure which for 
the parameters considered is a spiral galaxy-like one. 

In conclusion, we have shown that the initially ran- 
domly disturbed spinor condensate enters the equilibrium 
state in a step-like process. Counter intuitively, the dy- 
namics of the system becomes rich when it gets into the 
equilibrium. The atoms in components mp = +1 and 
mp = —1 spontaneously rotate in opposite directions 
and the magnetization shows the regular textures. More- 
over, the system exhibits coherence properties dependent 
on the strength of the spin-changing collisions. When 
such collisions are not strong enough but still allowing 
the transition to other Zeeman states only spinor com- 
ponents mp = ztl develop the mutual coherence. We 
observe that this coherence is periodically lost and re- 
stored. Simultaneously, the corresponding thermal parts 
get completely dephased. For larger values of coupling 
constants C2 the system enters a new regime where the 
coherence between components mp = and mp = ±1 
is developed and becomes as strong as for mp = +1 and 
mp = — 1 states. 
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